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The van der Waals interaction between two ground-state atoms is calculated for two electrically 
or magnetically polarizable particles embedded in a dispersive magneto-dielectric medium. Unlike 
previous calculations which infer the atom-atom interaction from the dilute-medium limit of the 
macroscopic, many-body van der Waals interaction, the interaction is calculated directly for the 
system of two atoms in a magneto-dielectric medium. Two approaches are presented, the first based 
on the quantized electromagnetic field in a dispersive medium without absorption and the second 
on Green functions that allow for absorption. We show that the correct van der Waals interactions 
are obtained regardless of whether absorption in the host medium is taken into account. 

PACS numbers: 12.20.-m, 42.50.Nn, 78.67.-n 



The van der Waals interaction between ground-state atoms in vacuum is often regarded as a consequence of the 
vacuum fluctuations of the electromagnetic field. It is well known that the interaction between two macroscopic, 
uncharged dielectric bodies cannot be obtained by pairwise addition of this interatomic van der Waals interaction 
except in the dilute-medium limit. In this limit Lifshitz, for instance, obtained the retarded and unrctarded pairwise 
van der Waals interaction between electrically polarizable atoms from a more general expression for the interaction 
energy of two dielectric half-spaces [l[ . The van der Waals interaction obtained in this way is identical to that obtained 
more directly by Casimir and Polder Q for the system of two atoms in vacuum. 

The continuing interest in Casimir effects and the related effects of vacuum field fluctuations in the case of dielectric 
media has stimulated some interest in the van der Waals interactions of atoms embedded in magneto-dielectric media. 
Recent work in this area takes essentially the original approach of Abrikosov et al. [|[ using Green functions 

for radiation in an absorbing medium and taking the dilute-medium limit to infer the atom-atom interaction. In 
this paper we obtain the atom-atom van der Waals interaction directly rather than from a dilute-medium limit of an 
interaction between macroscopic bodies, and we show that this interaction can be correctly obtained without explicit 
account of absorption in the host medium. 

We consider both electrically polarizable and magnetically polarizable atoms and present two derivations. The 
first, presented in Section [D] is based on the (electric or magnetic) dipole-dipole interaction induced in the two atoms 
by the "vacuum" field in the magneto-dielectric medium. For this purpose we employ simple expressions for the 
quantized electromagnetic field in a dispersive magneto-dielectric medium in which absorption is ignored 0]. The 
second, the subject of Section Hill employs Green functions and takes explicit account of absorption [n|. The fact 
that the same results are obtained regardless of whether absorption is accounted for appears to us to be of some 
interest, and physical reasons for it are given in Section llVl which also includes some further discussion and a brief 
summary of our results. 

II. VAN DER WAALS INTERACTIONS FROM QUANTIZED FIELD WITHOUT ABSORPTION 

We consider two identical atoms in a homogeneous and isotropic magneto-dielectric medium which we regard as 
a continuum with real electric permittivity e(uj) and magnetic permeability h(uj), and therefore real refractive index 
n(u). The calculation of the van der Waals interactions in this section will be based on the following expressions for 
the electric and magnetic fields in the non-absorbing magneto-dielectric (llj : 
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We employ a standard notation in which a,k\(t), a ] j cA (i) are Heisenberg-picture photon annihilation and creation 
operators for the plane-wave mode with wave vector k [|k| = k = n^cu^/c, = n(u>k)] and (linear) polarization unit 
vector ekA, [k-etA = 0, A = 1, 2]. The refractive index is = (e^/ife) 1 / 2 , where and Hk are the electric permittivity 
and magnetic permeability, respectively, at frequency u>k- The group index is 7fc = + uiudrik/ du>k, and V is the 
quantization volume for the box normalization of the plane-wave modes. The Hamiltonian is 

H = H A + H B + H F - d A (t) ■ E(r A , t) - d B (f) • E(r B , t) - m A (t) ■ H(r A , t) - m B (t) • H(r B , t), (3) 

where H A , H B , and Hp are the Hamiltonian operators for atom A, atom B, and the electromagnetic field, respectively. 
The electric dipole moment operators are denoted by d(i) and the magnetic dipole moment operators by m(i). 



A. Van der Waals interaction between electrically polarizable particles 



We first consider the van der Waals interaction between two electrically polarizable atoms separated by a distance 
R in the magneto-dielectric medium. The approach we will take follows closely that used by various authors for the 
derivation of the van der Waals interaction in the case in which the particles are in vacuum (n = 1) [l2|]. The basic 
idea of this method is that quantum vacuum fluctuations of the electromagnetic field in the medium induce in the 
atoms fluctuating electric and magnetic dipoles that interact with each other. The electric dipole moment induced in 
an atom at r by an electric field is d(t) = a e (£)E(r, t), where a e is the (real) electric polarizability of the atom. (Since 
it is only the real part of the polarizability that determines shifts in energy levels, we can assume without loss of 
generality throughout this paper that the polarizabilities of the guest atoms are real.) The dipole interaction energy 
between the two induced, fluctuating electric dipoles is then 



W ee (R) = £ ctf (w)af H (Ej (r A , kA)E, (r B , kA)> V% ( W> R), 
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where R = yb — v A , R = |R|. The two-point vacuum electric-field correlation function, summed over polarization 
states, follows easily from (flj: 



V(E l (r A ,kA)E,(r B ,kA)) = 2 ^Mfe (% _ , )e _i k .R 
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with ki = k,-/fc. The interaction potential V^ e (w, R) between two oscillating electric dipoles embedded in the magneto- 
dielectric is calculated as follows. From the Hamiltonian ^ and the Heisenberg equations of motion for the annihi- 
lation and creation operators, one obtains a formal expression for the electric field generated by the electric dipole B 
at the position of the electric dipole A: 
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where a = ds — (ds • R)R, b = d^ — 3(ds ■ R)R, and R = R/R. We have used the notation d(t) = dp(t), with d a 
unit vector specifying the direction of the electric dipole moment. We are interested in the case of two electric dipoles, 
both oscillating at frequency J: p B (t) = CBe^' 1 + C B e lJt and p A (t) = C A e~ lJt + C* A e lJt , where C A ,C B are 
arbitrary constants. We define the interaction V /ee (R) = — d A p A (t) ■ E e (r^, t). After performing the time integration 
for times t 3> and taking the time average of the resulting expression, we get 
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where as usual -q — > + . Here a = d A ■ ds — (d A ■ R)(d# • R) and b = d A ■ ds — 3(d A ■ R)(d_g • R). We can also define 
the interaction as V ee (i\) = — dsps(t) • E e (rs,t), which amounts to interchanging A and B above. This implies that 
we can take C a Cb — C* B C A . Note also that the time average p A (t)pB{t) = C A C* B + C A C B , so that we can write 
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Next we use the fact that n and (i (or actually their real parts which arc implicit here) are even functions of to to 
rewrite this as 
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Performing the trivial contour integrations, writing to instead of u>' for the dipolc frequencies, and using again k = 
n(u)u>/c, we obtain the electric dipole-dipole interaction tensor needed in Eq. (|4"|): 

V?f (w,R) = —-L[(Sij - 3RjRj)(cosfci? + kRsinkR) - (<% - RR,)/s 2 i? 2 cos kR}. (9) 



Using Eqs. (J3J) and ([5]), and passing to the continuum limit ^ k — > (V/8ir 3 ) dkk 2 J dQk = 
(V/8ir 3 c 3 ) J °° du>j LJ n 2 J uj 2 J dfl^, we obtain the van der Waals interaction energy between electrically polarizable par- 
ticles: 

Wee (-R) 



H 1 
^R 3 



uj 3 y?(uj) 
n{uj) 



kR sin 2kR + 2 cos 2kR, - 5 



sin2fci? cos2fci? 
6- 



sin2fci? 



kR 



k 2 R 2 



k 3 R 3 



(10) 

The integration path can be rotated using the fact that there are no there are no poles in the upper half of the complex 
plane. We obtain finally 
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where F(x) = x 4 + Ax 3 + 20x 2 + 48a; + 48. Recall that along the imaginary frequency axis the electric permittivity and 
refractive index are real and positive. The electric-electric van der Waals force that results from Eq. (fTTj) is therefore 
always attractive, regardless of the frequency dependence of e(w) and n(uj) . 

Our calculation based on the quantized fields (J]) and ^ gives a van der Waals interaction (fTTj) in full agreement 
with that obtained recently by Tomas for instance. The same is true for the other van der Waals interactions we 
calculate in this paper. The main point of this section is to show that correct results for van der Waals interactions 
involving ground-state atoms in dispersive media can be obtained straightforwardly, without having to go to a dilute- 
medium limit of an interaction between macroscopic bodies, and without having to introduce complexities arising 
from absorption. We discuss this further in Section IIV1 



B. Van der Waals interaction between magnetically polarizable particles 

We next use the same approach to calculate the van der Waals interaction between magnetic dipoles induced in 
the atoms by fluctuations of the zero-point magnetic field. For this the relation between an induced magnetic dipole 
moment at position r and the magnetic field is m(i) = a TO (t)H(r, t), where a m is the (real) magnetic polarizability 
of the atom. The dipolc interaction between the two induced, fluctuating magnetic dipoles is 

W mm (R) =^a^(o;)^( W )(H i (r A) kA)H i (r s ,kA))y^ m ( W ,R), (12) 

kA 

where the two-point vacuum magnetic field correlation function summed over polarizations is found from ([2]) to be 

^(H. t (r A ,kA)H,(r B ,kA)) = 2 * hn ^ {Sij _ k ikj ) e - ik R . (13) 

One can derive the magnetic dipole-dipole interaction tensor following steps similar to those above for the electric 
dipole-dipole interaction: 

T/" im (w,R) = — \[{5ij - 3R i R j ){coskR+ kRsmkR) ~ {8 l} - RR,)£; 2 i? 2 cos kR], (14) 
[ik R 

which differs from ([9]) simply by the replacement of €k by /ifc . The details of the evaluation of (p~2|) are essentially the 
same as for the electric van der Waals interaction and lead straightforwardly to the expression 
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Recall that along the imaginary frequency axis the magnetic permeability is real and positive. The magnetic-magnetic 
van der Waals force that results from Eq. (fT5)) is always attractive, regardless of the frequency dependence of fJ,(ui) 
and n(uj). 
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C. Van der Waals interaction between an electrically polarizable particle and a magnetically polarizable 

particle 

In calculating W ee (R) and W mm (R) it has not been necessary to account for the fact that the field operators in (j4|) 
and (|12p do not commute. Because of this noncommutativity, it is more appropriate to write W ee (R), for instance, 
in the symmetrized form 



Wee(R) = \Y, a * H a ? M [(E*(^,kA)Ej(r Bl kA)> + (E^, kA)E,(r A , kA))] V t f{u,TL) 

kA 

= Re^^MafH(E l (r A ,kA)E,(r B ,kA))^f( W ,R). (16) 

kA 

The forms (J5j and @, however, show that symmetrization is actually not required because the summation over k 
does not require us to distinguish between (E,(r^, kA)Ej(rg, kA)) and (Ej(rg, kA)Ej(r^, kA)). 

The situation in the case of the van der Waals interaction between an electrically polarizable particle and a mag- 
netically polarizable particle, however, is different because the electric-magnetic correlation function summed over 
polarization states, 

^<E,(r A ,kA)H J (r B ,kA)) = HZ^^g-** (17) 

A ^ k 

is not purely real when summed over k. (ey; is the Levi-Civita tensor.) Moreover the interaction tensor in this case, 
which we calculate to be 
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in a manner directly analogous to the electric-electric and magnetic-magnetic tensors, is antisymmetric. 

Let E m (r J 4, t) be the electric field (operator) at ya due to a magnetic dipole at yb- We write the interaction between 
the fluctuating electric and magnetic dipole moments in the symmetrized form 

W em (R) = -Rc^^( w )(E l (+) (r A ,kA)E r ( ni ) (r A ,kA)), (19) 

kA 

where ~E( + \rA,t) is the positive- frequency (photon annihilation) part of the source- free ("vacuum") electric field 
operator at r^, and eL ^(r^, t) is the negative-frequency (photon creation) part of the electric field produced by the 
magnetic dipole moment at yb- This electric field is induced by the source-free magnetic field H(re,£), so that the 
evaluation of (fT!)|) involves the electric-magnetic correlation function (jTTJ) . The calculation is essentially just the same 
as that presented by Farina et al [l3j for the case where the two particles are in free space, except of course that in 
our case the refractive index n(uj) appears: 



W em (R) = 4 ^ 4 / duu 2 af(iu)af n (iu)G 



2n(iu)uR 



e -2n(iu)uR/c ^0) 



where G(x) = (x + 2) 2 . n(iu) is real and positive, so that W em (R) is always repulsive, regardless of the frequency 
dependence of the refractive index. 



III. VAN DER WAALS INTERACTIONS FROM QUANTIZED FIELD WITH ABSORPTION 

In this section we will calculate the van der Waals interactions considered in the previous section for two atoms 
embedded in a magneto-dielectric, but now taking absorption in the host medium into account. We use the quantization 
procedure for the EM field in a dispersive and absorbing medium based on the Green- function formulation [l(i| . The 
dyadic Green function G(r,r',w) satisfies [l(| 



V x rc(r, w)V x j-e(r, u) 
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as well as the appropriate boundary conditions. Here k(t,uj) — ji 1 (r,ui). In an infinite, homogeneous material, 
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where R = r r', J? = |R|, R = R/i?, and k(uj) = n(u>)ui/c. The refractive index of the medium is given by 
n 2 (uj) = e((j)fx(u>), with e(u>) = e'(u>) + ie"(ui) the complex electric permittivity and = h'(uj) + ifj,"(ui) the 

complex magnetic permeability. 

The quantized electric field in dispersive, absorbing media may be written in the form (Iol | 

E(r,w)= J d 3 r'G x (r,r',Lu)-{ x (r',uj) + h.c., (23) 

A — e , m 

where the operators are bosonic operators satisfying the usual commutation relations: 

[f Xii (r,u ) ),ft, J (r',u>')}=6 XiX ,6 ij 6(u;-u;') ; [/ W (r,w),/ V)i (r'y)] = 0. (24) 

These operators may be regarded as being variables of the system composed of the EM field and the medium including 
the dissipative system. The electric and magnetic dyadic Green functions are defined in terms of the full Green function 

as 



G e (r,r» = *^/-Ime(r',^)G(r,r>), (25) 
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G m (r,r» = -i^/--Im K (r', W )[G(r 1 r' ) a))xV4 (26) 

Note that for an absorbing medium Ime(r, w) > 0, Im/x(r, u>) > 0, and Im/t(r, u>) < 0. The quantized magnetic field, 
similarly, may be written as 

H(r,w)= V /d 3 r'V r xG A (r,r',^).f A (r>) + /i.c., (27) 

and the total Hamiltonian for the free field is 

H = J J rfwfiwf !( r > w ) -h(r,uj). (28) 

A— e.rn ® 

We will require the following two-point vacuum field correlation functions obtained from these expressions: 

(E i (r I w)Ej(r',w , )> = " w'WG^y.w)]^, (29) 

J w c z 

(Hi(r,co)liUr',u')) = -n^^-^'WVrXV^xGfr.r',^, (30) 

(Ei(r )W )Hj(r>')> = -^ r T--%-^'MVr'xG(r,r', W )] iJ . (31) 

The electric and magnetic fields at the position of the atom A are given by the sum of the "vacuum" contributions, 
E (r j 4,w) and H (r^,w), plus the fields generated by the atom B which contain both electric dipole and magnetic 
dipolc components. An electric dipole located at position r# generates fields at position ya given by 

E e (r A ,w) = tu 2 af {uj)G(r A , r B ,uj) -Ho^b, lu) + h.c, (32) 
H e (r j4 ,w) = -iujcn(uj)af(uj){\7 rA x G{y a ,y b ,uj)} ■ E (r B ,w) + h.c, (33) 

while the fields generated by a magnetic dipole at position yb are 

E m (r A ,w) = -iwrt(w)a^(w)[V rA x G(r A ,r B ,u)] -Ho^u;) +/i.c, (34) 
H m (rA,w) = cK(w)fc 2 (w)a^(w)G(r j4 ,r B ,w) • H (r B ,w) + /i.e. (35) 
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The vacuum expectation value of the van der Waals energy may be written as the sum of three contributions, 
one purely electric, one purely magnetic, and one mixed. Using the above expressions for the two-point correlation 
functions of the EM field, one can easily find each of these terms. The purely electric part stems from the p • E 
interaction, and is found to be 
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Comparing this expression with Eq.Q we see that o-> 2 Im[G(ivi, r_B, ^>)]ij is related to the electric dipole-dipole interac- 
tion tensor U^ e (u;, R), and that u> 2 Re[G(r A , r^, u>)]ij is related to the solid-angle integration of the two-point vacuum 
electric-field correlation function summed over polarization states, given in Eq.((5|). After rotation in the complex 
plane (u> — > iu), we can re- write this expression as 

W ee (R) = -— duai(iu)af(iu)u 4 Tr[G{r A ,r B ,iu)-G(r A ,r B ,iu)}. (37) 
27T Jo 

The purely magnetic part, similarly, comes from the m • H interaction: 
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where Eq. (|21|l and a rotation in the complex plane were used in obtaining the last equality Finally, there are two 
electric-magnetic terms, one arising from the p • E interaction, and one from the m • H interaction. They result in 
the mixed interaction 
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The traces appearing in the integrands of Eqs. (|37l38l39p can be explicitly computed given the form of the dyadic Green 
function evaluated at the imaginary frequency w = iu. The final result for the complete van der Waals interaction 
energy between two ground state atoms embedded in an absorbing and dispersive medium is then 
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where again F(z) — z 4 + 4z 3 + 2Qz 2 + 48z + 48 and G(z) = (z + 2) 2 . This is identical to the complete van der Waals 
interaction obtained in Section HT1 



IV. DISCUSSION 



Since the van der Waals interaction between electrically polarizablc particles is the most important, a rough model 
for the modification of the vacuum interaction by the medium might be of interest. Let us consider a two- level model 
in which the polarizabilities are 

aft W ) = a?( W ) S a( W ) = 53^, (41) 
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FIG. 1: The van der Waals interaction (|lip divided by the London interaction for a two-level model. The upper curve is for 
two atoms in vacuum (C = 0), the lower curve for two atoms in a dielectric defined by (|45|) with C — 3. 



where d and uj$ are respectively the (real) transition electric dipole moment and the transition angular frequency, and 
similarly 

n{uo) = [1 + 4:nNa(uj)} 1/2 , (42) 

where N is the atomic density of the host medium and we take \i = 1. It is useful to normalize (jllj) to the familiar, 
nonretarded London form of the interaction: 

We define 
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where r = ujqR/c and e and n are evaluated at iuJoy: 
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e(iu y) = n 2 (iuj y) = 1 + (45) 
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where C = 8irNd 2 /SHujq. In the limit C = and r — ► 0, D — ► 1; for j" — ► oo, £) — > 23/37rr, or W ee = -23ftca 2 (0)/47ri? 7 , 
the famous Casimir-Polder result. Figure 1 plots D(R) for C = (vacuum) and C = 3. 

This simple model is not in any sense meant to be a realistic, quantitative description of the van der Waals 
interaction between two atoms embedded in a dielectric. It does, however, suggest that the predominant effect of the 
host medium on the van der Waals interaction is to weaken it, without substantially changing the distance dependence 
in either the nonretarded or retarded regimes. More realistic models of the van der Waals interaction in a liquid, for 
example, must take into account local field corrections, as has been discussed, for instance, by Abrikosov, et al. Q 



and McLachlan [14 1 



Our results for the electric-electric, magnetic-magnetic, and electric-magnetic van der Waals interactions between 
two atoms embedded in a dispersive magneto-dielectric medium are in agreement with those obtained previously 
3"0- Unlike previous derivations, however, we obtained the interactions directly rather than inferring them from 
the dilute-medium limit of van der Waals interactions between macroscopic bodies, where the interactions can be 
obtained by pairwise summations of interatomic interactions [H-Q. 

What is perhaps more interesting, however, concerns the role of absorption in the host medium. In the preceding 
section we accounted for dissipation (absorption) in the host magneto-dielectric, as in previous work. In Section HT1 
however, we ignored any possibility that the medium could be absorbing: we worked with expressions for the electric 
and magnetic fields that derive directly from the assumption that the medium is non-absorbing [llj . 

The same situation holds, for instance, in the case of the Lifshitz formula for the van der Waals interaction between 
infinite, plane-parallel dielectric media. Lifshitz's original derivation, and various derivations that followed, include 
the imaginary (dissipative) part of the dielectric function e(w). Derivations of the Lifshitz formula based on changes 
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in zero-point field energy arising from the dielectric media, however, make no reference to dissipation [15j |. [161 ]. [171 ]. 
In other words, the Lifshitz formula can be derived without explicit accounting for absorption. 

Ginzburg [Tj[ has also noted that (macroscopic) van der Waals interactions can be correctly derived based on changes 
in zero-point field energy, without accounting for absorption. He remarks, in connection with such derivations [l6| . 
that "oddly enough there is no mention that they consider directly only transparent media," and then gives reasons 
why the van der Waals (electric-electric) interaction for real media can be obtained by presuming non-absorbing 
media: "Firstly, the permittivities ... are functions. Secondly, the function e(w) is always real on the imaginary axis." 
The results (|11[) , (|15p . and (|20p . for instance, all involve permittivities and permeabilities on the imaginary axis. 

A simple and more physical explanation can be given for why (ground-state) van der Waals interactions calculated 
for non-absorbing media apply directly to real (absorbing) media, as we have found for the electric-electric, magnetic- 
magnetic, and electric-magnetic van der Waals interactions between two atoms embedded in a magnetic-dielectric 
medium. At zero-temperature, for instance, any atom of the host medium is in its ground state and can absorb 
radiation that is resonant with one of its transitions to an excited state. It cannot, of course, absorb from the vacuum 
field: in this case the fluctuations in the field that might induce absorption are exactly cancelled by fluctuations in 
the atom itself [13]. The same is true when the atom is part of a magneto-dielectric in which is embedded, as in 
examples considered in this paper, two guest atoms. Any atom of the host medium still finds itself in a vacuum field 
state, regardless of the nature or the number of guest atoms. The host and guest atoms modify the modes of the 
field from the simple plane waves of a pure vacuum, but the field remains in a vacuum state |vac) (f\^\v&c) = for 
any mode). Just as in free space, therefore, there is no absorption unless it is possible to populate one or more field 
modes; this would be the case only if there were an applied external field or if one or more atoms is excited, leading 
to the possibility that a different atom could absorb its emitted (real) photon. In other words, absorption by the host 
medium would play a role if we were to consider a van der Waals interaction involving excited atoms. Otherwise 
one can expect to obtain correct van der Waals interactions without having to account for the absorption that is 
always present in a real medium. This expectation applies, of course, regardless of how many atoms are involved and 
regardless of the shape of any macroscopic bodies for which the van der Waals forces are to be calculated. 

In light of recent interest in negative- index media [il^ . [20| , let us reconsider specifically the most important of the 
van der Waals interactions we have calculated in this paper, namely that between two electrically polarizable atoms. 
First we note that the expressions for the quantized fields in Section [IT] are directly applicable to negative- index media 
[Tl| . In a negative index medium n, e, and fi are all negative at some frequency or range of frequencies. This would 
at first glance suggest that the dipole-dipole interaction V^ e (u),TC) [Eq. ©] changes sign at frequencies for which the 
refractive index is negative. To see that this is not the case, note that Eq. © is unchanged whenever n(uj), e(u>), and 
all change sign within any frequency range. This means that the dipole-dipole interaction in a negative-index 
medium does not change sign, and in particular that is directly applicable in general provided we just replace 
n(ui), e(u>), and fi(u>) by their absolute values. It follows similarly that the van der Waals interaction does not change 
sign or undergo any other significant change in a magneto-dielectric medium. The same conclusion applies to the 
magnetic- magnetic and mixed van der Waals interactions, and is in agreement with the conclusions of Buhmann et 
al. [|. 

The pairwise electric-electric and magnetic-magnetic van der Waals interactions are always attractive, whereas the 
pairwisc electric-magnetic interaction is always repulsive. These results apply also in the case of negative-index media, 
at least to the extent that such media can be modeled as continua. It is well known, however, that non-pairwise van 
der Waals interactions can be repulsive [2l[. Evidently repulsive Casimir effects such as those recently suggested by 
Henkel and Joulain [22| and by Leonhardt and Philbin [23| in the case of negative-index media must in some way 
involve either non-pairwise interactions or electric-magnetic van der Waals interactions. 
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